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Abstract 

We give a fully polynomial randomized approximation scheme to 
compute a lower bound for the matching polynomial of any weighted 
graph at a positive argument. For the matching polynomial of complete 
bipartite graphs with bounded weights these lower bounds are asymp- 
totically optimal. 
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1 Introduction 

Let G = (V, E) be an undirected graph, (with no self-loops), on the set of 
vertices V and the set of edges E. A set of edges M <Z E is called a matching 
if no two distinct edges ei,e2 G M have a common vertex. M is called a 
k-matching if j^M = k. For k G N let Ai^G) be the set of /c-matchings in G. 
(Aik(G) = for k > |_4pJ.) If #V = 2n is even then an n-matching is called 
a perfect matching. 

Let uj : E — > (0, oo) be a weight function, which associate with edge e G E 
a positive weight u(e). We call G^ = (V, E, u) a weighted graph. Denote by i 
the weight l : E — ► {1}. Then G can be identified with G^. 

Let M G M.k{G). Then the weight of the matching is denned as uj(M) : = 
Y\ e< - M ui(e). The total weighted k-matching of is defined: 

<j>{k,GJ) := u(M),keN 

MeM k (G) 
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where (j){k,GJ) = if Aik{G) = for any k e N. Furthermore we let 
0(0, G u ) := 1. Note that 4>(k, G b ) = #M.k(G), i.e. the number of fc-matchings 
in G for any k G N. The weighted matching polynomial of G^ is defined by: 

G w ) := £ n = L^J • 

fc=0 

This polynomial is fundamental in the monomer-dimer model in statistical 
physics [3l [12], and for u = 1 in combinatorics. Note that if #y is even 
then $(0, G u ) is the total weighted perfect matching of G. (Some authors 
consider the polynomial tL^-J$(t^ 1 ) QJ} instead of $(t, G u ).) It is known that 
nonzero roots of a weighted matching polynomial of G are real and negative 
[T2] . Observe that $(1, G L ) the total number monomer-dimer coverings of G. 

Let G be a bipartite graph, i.e., V — V\ U V 2 and E 1 C V\ x V2. In the 
special case of a bipartite graph where n = $V\ = #V 2 , it is well known that 
4>(n,G) is given as perm B(G), the permanent of the incidence matrix B(G) 
of the bipartite graph G. It was shown by Valiant that the computation of 
the permanent of a (0,1) matrix is #P-complete [17]. Hence, it is believed 
that the computation of the number of perfect matching in a general bipartite 
graph satisfying #Vi = jfV 2 cannot be polynomial. 

In a recent paper Jerrum, Sinclair and Vigoda gave a fully-polynomial ran- 
domized approximation scheme (fpras) to compute the permanent of a non- 
negative matrix [13]. (See also Barvinok [1] for computing the permanents 
within a simply exponential factor, and Friedland, Rider and Zeitouni [pj for 
concentration of permanent estimators for certain large positive matrices.) [13] 
yields the existence a fpras to compute the total weighted perfect matching in 
a general bipartite graph satisfying #V\ = j^V 2 - In a recent paper of Levy and 
the author it was shown that there exists fpras to compute the total weighted 
fc-matchings for any bipartite graph G and any integer k £ [1, ^]. In particu- 
lar, the generating matching polynomial of any bipartite graph G has a fpras. 
This observation can be used to find a fast computable approximation to the 
pressure function, as discussed in [8], for certain families of infinite graphs 
appearing in many models of statistical mechanics, like the integer lattice Z d . 

The MCMC, (Monte Carlo Markov Chain), algorithm for computing the 
total weighted perfect matching in a general bipartite graph satisfying $V\ = 
#V 2 , outlined in [13], can be applied to estimate the total weighted perfect 
matchings in a weighted non-bipartite graph with even number of vertices. 
However the proof in [13], that shows this algorithm is frpas for bipartite 
graphs, fails for non-bipartite graphs. Similarly, the proof of concentration 
results given in [9] do not seem to work for non-bipartite graphs. The tech- 
nique introduced by Barvinok in [T] to estimate the number of weighted perfect 
matching in bipartite graphs, does extend to the estimate of total weighted 
perfect matchings in a general non-bipartite graph with even number of ver- 
tices, when one uses real or complex Gaussian distribution. (See the discussion 
in §5.) 
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In this paper we give a fpras for computing a lower bound $(£, GJ) for the 
weighted generated function <3>(t, GJ) for a fixed t > 0. We show that this 
lower bound has a multiplicative error at most exp(jV"min(|r, Ci)), see (ll.7p . 
where a 2 is the maximal weight of edges of G and 

If 2 

Ci = = / \og{x 2 )e~^dx = 1.270362845 . . . (1.1) 

V27T 

These estimates are similar in nature to heuristic computations of Baxter [2], 
where he showed that his computation for the dimers on Z 2 lattice are very 
precise away from only dimer configurations, i.e. perfect matchings. (The 
results of heuristic computations of Baxter were recently confirmed in [§].) 
We show that that for the matching polynomial of complete bipartite graphs 
with weights in [b 2 , a 2 }, < b < a, this lower bound is asymptotically optimal. 

We now describe briefly our technical results. With each weighted graph 
G w associate a skew symmetric matrix A = [ay]^. =1 G R NxN , A T = —A, 
where iV := jJV , as follows. Identify E with (N) := {1,...,N}, and each 
edge e G E with the corresponding unordered pair ^ j G (N). Then 

dij if and only G E. Furthermore for 1 < i < j < N, G E 

dij = For 1 < i < j < N let Xij be a set of (^) independent 

random variables with 

E = 0, E x\ = 1, 1 < i < j < N. (1.2) 

Let x := (xn, . . . , xuv, ^22, • • • , x nn)- We view x as a random vector vari- 

/JV+1\ 

able with values £ = (£n, . . . , £,nn) £^ 2 Let be the following skew- 
symmetric random matrix 

Ya '■— [Oij^min(ij') max(«,j)]j j'=l G M. . (1-3) 

A variation of the Godsil-Gutman estimator [TU] states 

E det( v / t/v + Y A )) = G w ) if = #"1/ is even, (1.4) 
E det(VtI N + Y A ) = Vi<S>(t, GJ) li N = #V is odd. (1.5) 

for any t > 0. Here Zjy stands for N x N identity matrix. 
We show the concentration of logde^v^Tv + Y A ) around 

log$(t,G w ) :=E logdet( v / t/7v + Y A ) (1.6) 

using [TT]. These concentration results show that $(t, GJ) has a fpras. Jensen 
inequalities yield that $(t, (j w ) < $(t, G w ). Together with an upper estimate 
we have the following bounds: 

llog<§(t,Gy < l]og$(t,G w ) < llog$(t,G w ) + min(^,C 1 ) (1.7) 

where a = max |ay|. The above inequality hold also for t = 0. (For A^ even 
and t = this result is due to Barvinok pQ §7].) It is our hope that by refining 
the techniques we are using one can show that $(£, GJ) has fpras for any t > 0. 
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2 Preliminary results 



Lemma 2.1 Let G = (V, E) be an undirected graph on N vertices. Let 
uj : V —* (0, oo) be a given weight function. Let A = —A T G M. nxn be the 
corresponding real skew symmetric matrix defined in §1. Assume that Xij,i = 
1, . . . ,j,j = 1, . . . , N are ( Ar ^ _1 ) independent random variables, normalized by 
the conditions U.!H\) . Let Ya G M. NxN be the skew symmetric real matrix defined 
by (CLf. Then (I-MP-Z.5]) hold. 

Proof. Let \fi = s. Observe first that det(siAr + Ya) is a sum of AT! 
monomials, where each monomial is of degree at most 2 in the variables Xij 
for i < j and of degree m invariable s. The total degree of each monomial is 
N. The expected value of such a monomial is zero if at least the degree of one 
of the variables x^ is one. So it is left to consider the expected value of all 
monomials, where the degree if each x^ is or 2, which are called nontrivial 
monomials. 

Assume first that A" is even. Observe that if a monomial contains s of 
odd power than it must be linear at least in one x^. Hence its expected value 
is zero. Thus E det (sijv + Ya) is a polynomial in s 2 . Consider a nontrivial 
monomial such that the power of s is A" — 2m. Note that this monomial is of 
the form rs N ~ 2m Yl(i j)^M j)) x %i f° r some m matching M G Ai m . Here 
(— l) m r is the sign of the corresponding permutation o : (N) — > (N). Since 
a(i) = j, a (J) = i for any edge (i,j) G M, and a(i) = i for all vertices i which 
are not covered by M we deduce that r = 1. Hence the expected value of this 
monomial is s N ~ 2m Y[ e& M u ( e ) ■ This proves (11 Ah . The identity ( 11. 5ft is shown 
similarly. □ 

Recall the following well known result: 

Lemma 2.2 Let A = —A T G R 7Vx7V be a skew symmetric matrix. Then 
B := lA, where l := y— 1, is a hermitian matrix. Arrange the eigenvalues of 
B in a decreasing order: \\{B) > . . . > Xn(B). Then 

X N - i+1 {B) = -\(B) for i = 1, . . . , N. (2.1) 

In particular 

N 

det{Vil N + A) = Y[ y/t + Xi(B) 2 . (2.2) 
t=i 

Proof. Clearly, B is hermitian. Hence all the eigenvalues of B are real. 
Arrange these eigenvalues in a decreasing order. So —i\j(B),j = 1,...,N 
are the eigenvalues of A. Since A is real valued, the nonzero eigenvalues of A 
must be in conjugate pairs. Hence equality (12. ip holds. Observe next that if 
X k (A) = -iX k {B) ^ then 

(Vt + X k (A))(Vt + X N - k+1 (A)) = ^t + X k {B)^t + X N - k+l {By. 
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As the eigenvalues of \/il N + A are \ft + Xk(A), k = 1, . . . , N we deduce (12. 2p . 

□ 



3 Concentration for Gaussian entries 

In this section we assume that each normalized real Gaussian variable, 

i.e satisfying (11.21) . Recall that a function / : R — > R is called Lipschitz 
function, or Lipschtzian, if there exists L G [0, oo) such that ^^-^[^ < L for 
all x 7^ y G R. The smallest possible L for a Lipschitz function is denoted by 
\f\ c . Let A N C R nxn , iA N C C nxn denote the set of Ax A" real skew symmetric 
matrices, and the set of N x N hermitian matrices of the form iA, A G An- 
With each A G A N we associate a weighted graph = (V,E,uj), where 
V = (A), G V aij ^ 0, uj((i,j)) = \aij\ 2 . Denote by a := max|ajj|. 

To avoid the trivialities we assume that a > 0. Note that a 2 is the maximal 
weight of the edges in G u . Let Ya be the random skew symmetric matrix given 
by (11.31) and denote by Xa the random hermitian matrix Xa '■= ^p 1 ^- 

Let / : R — > R be a Lipschitz function. As in [TTJ consider the following 
F : iAn — > R given by the trace formula: 

1 - 

F{B) = trjy f(B) := - ]T 5 G iA N . 

Denote by E tr N (f(XA)) the expected value of the function tr A r(/(A J 4)). The 
concentration result [TT1 Thm 1.1(b)] states: 

N 2 r 2 

Pr(| tr^(/(A A )) - E ti N (f(X A ))\ > r) < 2e _i ^I (3.1) 

(Recall that the normalized Gaussian distribution satisfies the log Sobolev 
inequality with c = 1.) We remark that since the entries of Xa are either zero 
or pure imaginary one can replace the factor 8 in the inequality (13.11) by the 
factor 2. See for example the results in [151 8.5]. 

Lemma 3.1 Let ^ A = [ay] G An, a = max|aij|,t G (0,oo) ; Xy, 1 <i< 
j < N be independent Gaussian satisfying U.S\) . Let Ya G An be the random 
skew symmetric matrix given by U.3\) . Then 

Pr(| logdet( v / t/jv + Y A ) — E logdet( v / t/iv + Y A )\ > Nr) < 2e~^ . (3.2) 
Proof. Let ft(x) := |log(-^ + x 2 ). f t is differentiable and 

\(ft)c\ — max I ft( x )\ — — F- 
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Apply (13. ip to f t . Observe that the right-hand side of (13.11) is equal to the 
right-hand side of (13. 2p . Use (12. 2p to deduce that 



N /— N 



N tr s{ft (X A ) ) = £ log J 1 + A, {X A f = J2 log J ^ + ' ' ' * 1 1 * 



, N ^— ' \ N N 

i=i i=i 



1 " 1 

-N log TV + log J] \/tTpM = --iV log N + log det(V*Ijv + Y A ) . 



i=l 



Hence the left-had sides of (13.11) and ( 13. 2ft are equivalent. □ 

The following lemma is well known, e.g. [91 p'1566], and we bring its proof 
for completeness. 

Lemma 3.2 Let U be a real random variable with a finite expected value 
E U. Then e EU < E e u . Assume that the following condition hold 

Pi(U - E U > r) < 2e- Kr2 for each r e (0, oo) and some K > 0. (3.3) 

Then 

e EU <Ee u < e EU (l + ^L). (3.4) 

Proof. Since e u is convex, the inequality e EU < E e u follows from 
Jensen inequality. Let fi := E U and F(u) := Pr(U < u) be the cumulative 
distribution function of U. We claim that 

Ee y <e"+ / e u (l - F{u))du. (3.5) 

J fl<U 

Clearly 



Ee u = I e u dF(u)= I e u dF(u) + j e u dF(u). (3.6) 



1 U<fl J fJ,<U 

Since e u < e^ for u < /i we deduce that 

I e u dF{u) < ef*F(fi). 

J U<fl 

We now estimate the second integral in the right-hand side of ( 13.61) . Recall 
that F(u) is an nondecreasing function continuous from the right satisfying 
F(+oo) = 1. Hence e u (F(u) - 1) < for all u G R. For any R > fi use 
integration by parts to deduce 

/ e u dF(u) = e u {F(u) - 1)|J+ / e u (l - < 

J H<u<R J \i<u<R 

e"(l-F(^))+ / e s (l -F(«))du. 



So 



e u dF(u) < e M (l - F(fi)) + / e u (l - F(u))du, 

and (13.51) holds. 

Assume now that ( 13. 3D holds. Thus 

1 - F(u) = Pr(C/ > u) < 2e^ { "^ )2 for any u > fi. 

Hence 

/ e u (l - F{u))du <2 I e u - K{u -» )2 du < 

J fl<U J fl<U 

2e» r e- K ^-^^du 2e " e " 



VK 



TV 



Combine the above inequality with (13.51) to deduce the right-hand side of (13. 4p . 

□ 



Corollary 3.3 Let the assumptions of Lemma \3. 1\ hold. Then 

1 1 1 1 V%Nae^ 

■ log < -log < -log$(t,G u ) + -log(l + 



4 FPRAS for computing log$(t, 

Let £ G M 7Vx7V . For jfe G N denote by ® k B G M fcAfxfcAr the block diagonal 
matrix diag(5, . . . ,B). {®kB is a direct sum of k copies of B.) Note that if 

k 

B G An then © fc S G A kN . Clearly, 

det(s47v + ®kB) = (det(s/7v + B)) k for any B G M JVxAr and s G E. (4.1) 

Let A G An, and Ya be the random matrix defined by (11.31) . By Ya(£) we 
mean the skew symmetric matrix [ajj£ min (jj) max (jj)]^ =1 , which is a sampling of 
Fa- Let Xij,l < i < j < kN be ( fciV 2 +1 ) normal G aussian independent random 
variables. Consider the random matrix Y® h A- Then a sampling 

2 } =diag(y A (^),...,F A (^)),e t GR( =» ),i = 
is equivalent to k sampling of Ya- 

Theorem 4.1 Let ^ A = [a^] G An, a = max|aij|,t G (0,oo) ; x^-, 1 < 
i < j < N be independent Gaussian satisfying $1.21) . Let Ya G An be the 
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random skew symmetric matrix given by U.3\) . Let Y^^), . . . , Y A (£ h ) be k 
samplings ofY A . Then 

i k 2 

I r — ■» — ~ tkNr 

PT(\-J2^gdet(Vil N + Y A (^)) -log$(t,G u )| > Nr) < 2eSS^. (4.2) 

i=i 

In particular the inequality 

i tog < ilog$(t,G w ) < ilog$(t,G u ) + ^ (4.3) 

Hence an approximation of&(t,G w ) by (J^iLi det(\/tijv + 
fully-polynomial randomized approximation scheme. 

Proof. Use ( fl~T]) to obtain 

bgdet(V*4iv + *© fc A(£)) = ^logdet( v / t/iv + Y A (Q) 

i=i 

Hence 

E logdet( v / t4iv + Y @hA ) = kE \ogdet((Vil N + Y A ) = Hog<£(t, G w ) (4.4) 
Apply (13. 2p to Y^ kA to deduce (14. 2p . Observe next that 

E det(VtI k N + Y^a) = E det((VtI N + Y A f = $(t, G w f. (4.5) 
Use Lemma [3.21 for the random variable log det(\/tIkN + Y® kA ) to deduce 

llog$(t,G w ) < lbg$(t,G w ) < llog$(*,G u ) + 



1 , s/SkNae' 
+— bg(l 



Let A; — > oo to deduce (14. 3p . 

We now show that (14.21) gives fpras for computing $(£, G w ) in sense of 
Let e, 5 G (0,1). Choose 

r- e fc -r^Mi 

r "2iV' ^ 

Then 

Pr(l - e < d6t( ^ + Y ^ < 1 + e) > 1 - i. 
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Observe next that 



Prflx^^log— )<— . 

Hence with probability 1 — | at least, the absolute of each off-diagonal of 
i = 1) ■ ■ ■ j k is bounded by a-v/ 2 log In this case all the entries of 

y/il^ + yk(£i)) ar e polynomial in a, a/^, AT, -, log The length of the storage 
of each entry is logarithmic in the above quantities. 

Finally observe that we need 0(N 3 ) to compute de^-v/tiTv + ^A Hence 
the total number of computations for our estimate is of order 

rVivV 2 log r 1 . 

□ 

The quantity ilog<I>(£, G u ) can be viewed as the exponential growth of 
log$(i, G u ) in terms of the number of vertices N of G. Note that since the 
total number of matching of a graph G is given by $(l,G t ), Theorem 14.11 
combined with (jl.7p yields that the exponential growth of the computable 
lower bound $(1,GJ differs by | at most from the exponential growth of 
$(1, G L ). Note that for complete graphs on 2n, the exponential growth of the 
number of perfect matching matchings is of order log2n — 1. For fc-regular 
bipartite graphs on 2n vertices the results of [H [7] imply the inequality that 
for n big enough the exponential growth of the total number of matchings 
is at least log A; — 1. Thus for graphs G on 2n vertices containing, bipartite 
fc-regular graphs on 2n vertices, with k > 5 and n big enough, $(1, G t ) has a 
positive exponential growth. 

5 Another estimate of log$(t, GJ) — log$(t, GJ) 

Lemma 5.1 Let X be a real Gaussian random variable. Then 

log E X 2 — E log X 2 < Ci, (5.1) 

where C\ is given by M.l\) . Equality holds if and only i/EI = fl. 

Proof. Clearly, it is enough to prove the lemma in the case X = Y + a, 
where Y is a normalized by fll.2p and a > 0. In that case the left-hand side of 
(15. ip is equal to 

If x 2 

g(a) := log(l + a 2 ) -==. \ log((x + a) 2 )e~~ dx. 

V27T Jr 

We used the software Maple to show that f(a) is a decreasing function on 
[0, oo). So /(0) = C\ and lim^oo f(a) = 0. This proves the inequality (15. 11) . 
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Equality holds if and only if X = bY for some b 7^ 0. 



□ 



Denote by S n C IR nxn the space of n x n real symmetric matrices. A 
polynomial P : R n — > R is of degree 2 if 

P(x) = x T Qx + 2a T x + 6, 
x = (x u . . . , x n ) T , a = (ai, . . . , a n ) T e R n , Q e S n , b e R. 

(We allow here the case (5 = 0.) The quadratic form P^ : R n+1 — > R induced 
by P is given 

P h {y) = y T Q h y,Q h = 

Clearly, P(x) = P/ l ((x T , 1) T ). P is called a nonnegative polynomial if P(x) > 
for all x G R n . It is well known and a straightforward fact that P is non- 
negative if and only if Qh is a nonnegative definite matrix. 

The following lemma is a generalization of [U Thm 4.2, (1)]. 

Lemma 5.2 Let P : R n — > R fre a nonzero nonnegative quadratic poly- 
nomial. Let X\, . . . , X n be n- Gaussian random variables, and denote X : = 
(X u ...,X n ) T . Then 

E log P(X) < log E P(X) < E log P(X) + Ci, (5.2) 

where C\ is given by U.l\) . 

Proof. We may assume without a loss of generality that E P = 1. In 
view of the concavity of log we need to show the right-hand side of (15.21) . Since 
Qh is nonnegative definite it follows that 

m 

P(x) = M a 7* + bi) 2 , a; G R n , h G R, A, > 0, i = 1, . . . , m, 

m 

E (a t T X + 6j) 2 = 1, z = 1, . . . m, ^A i = l. 

i=l 

Note that one can have at most one a^ = 0, and in that case then bf = 1. The 
concavity of log yields 

TO 

logP(X)>^A i log(a7X + 6 i ) 2 . 

i=l 

(We assume that logO = —00.) Note that if aj 7^ then a,X + b{ is Gaussian. 
Lemma O yields E logP(X) > -C 1 . □ 



Q a 

a T b 



G S n+ i,y = (yi, . . . ,y n +i) 
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Theorem 5.3 Let the assumptions of Theorem 4-1 hold. Then ( 1.1) holds. 
Proof. In view of (14.31) it is left to show 

log $(£, G u ) < log G u ) + (N- 1)C X . (5.3) 



Let A = [ajj]^ =1 G A N . Recall that detA = (pfaf A) 2 , where pfaf A is the 
pfaffian. (So pfaf A = if n is odd.) Let a^ = (a-n, . . . , a(j_i)j) T G 1R* _1 , z = 
2, . . . , n. We view pfaf A as multilinear polynomial Pf (a2, . . . , a n ) of total de- 
gree |, which is linear in each vector variable a«. (Any polynomial of noninteger 
total degree is zero polynomial by definition.) 

Denote by Qk, n the set of subsets of (n) of cardinality k G [l,n]. Each 
a G Qk, n is viewed as a — {ix, . . . 1 < h < ■ ■ ■ < ik < rn. For any matrix 
B = [bij] G R nxn and a G Q k , n we define G R kxk as the principal 

submatrix [b ai aj]i t j=v Then for A = [a^] G A n denote 

Pf a (a 2 , . . . , a n ) := pfaf A[a|a]. 

Then Pf Q (a 2 , . . . , a n ) is a multilinear polynomial of total degree |, which is 
linear in each a,. Hence 

det(sI N + A) = s N + J2 sN ' k Yl Pfa(a 2 ,...,a w ) 2 , for any A G A N . (5.4) 

k=l «eQ fci „ 

View aj G W~ l as a variable while all other a 2 , . . . , a^r are fixed. Then for 
s > the above polynomial is quadratic and nonnegative. Group the (^) 
independent normalized random Gaussian variables A^, 1 < i < j < N into 
N — 1 random vectors Xj := (Xu, . . . , X^_i^) T , i = 2, . . . , N. Consider now 
Y A . Let 

P{X 2 ,...,X N ) :=det{VtI N + Y A ) t>0. 

Then P(X 2 , . . . ,Xn) is a nonnegative quadratic polynomial in each Xj,j = 

2, . . . , N. Denote by E j the expectation with respect to the variables Xu, . . . , Xu^iy. 

(E3D yields that 

Fj(X 2 , . . . , Xj) := E j +1 . . . E N F(X 2 , . . . , X^r) 

is a nonnegative quadratic polynomial in each Xj, j = 2, . . . , i. Lemma 15.21 
yields 

logE iPi(X. 2 , . . . ,Xi) < E ilogFj(X 2 , . . . ,X<) + Ci, % = 2, . . . , N. 
Hence 

log$(t, G u ) = logE 2 P 2 (X 2 ) < E 2 logP 2 (X 2 ) + d < 
E 2 E 3 logP3(X 2 ,X 3 )+2C 1 <...< 
E 2 E 3 . . . E N log P(X 2 , X 3 , . . . , Xjv) + (AT - l)d = 

log<£(i,Gy + (A-1)C' 1 . 

□ 
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6 Bipartite graphs 



Assume that G = (V, E) is a bipartite graph. So V — V1UV2, E C E\ X E 2 and 
N = m + n. Assume for convenience of notation that m : j^V\ < n := #V*j. 
Thus E C (m) x (n), so each e G E is identified uniquely with (i, j) G (m) x (n). 
Let C = [cy] G M mxn be the weight matrix associated with the weights u> : 
E — > (0,oo). So c.y = if G" E and c^- = y/oj{i, j) if (i,jf) G £7. Let 
Xij,i = 1, ... ,m, j = 1, ... ,ra be ran independent normalized real Gaussian 
variables. Let Vc ='■ [cij%ij] £ ]R mx?1 be a random matrix. Then the skew 
symmetric matrix A associated with G^ is given by and the corresponding 
random matrices Ya, Xa are given as 



A = 
Denote by 



C 
-C T 



Ya 



U c 

-ul 



■ X A = Y A . (6.i; 

\Jm + n 



(Tx{V c )>...>a m {V c )>0 (6.2) 

be the first m singular values of Vc- Then the eigenvalues of Y A consists of 
n — ra zero eigenvalues and the following 2m eigenvalues: 

±ia 1 {U c ),...,±uT m (U c ). (6.3) 

Hence 

m 

det(Vil m+n + Y A ) = J](t + a t (V c ) 2 ). (6.4) 

8=1 

In [9] the authors considered the random matrix Vc '■= VcV^ G IR mxm . Note 
that the eigenvalues of Vc are 

vI(V c )>...>a 2 m (V c ). (6.5) 

Furthermore, one has the equality E detVc = 4>(m,G w ). Let K m<l% be the 
complete bipartite graph on V\ = (m) , V 2 = (n) vertices. Assume that 1 < 
m < n. Let < b < a be fixed. Denote by f2 m n) [ 6 2 )a 2] the sets of all weights 
lu : (m) x (n) — > [6 2 ,a 2 ]. Recall that each u G £l m , n ,[& 2 ,a 2 ] induces the positive 
matrix C(oj) = [cij(u)] G R mxri , where Cij(co) G [b, a]. It was shown in [9] that 
- log det Vc(cj) concentrates at - log (f>(m, K m>n ^) with probability 1 as n — ► 00. 
More precisely 

limsup sup Pr( — I log det Vc{j) - log0(ra, if m , n ,w)| > 5) = (6.6) 



,[62, a 2] n 



for any 5 > 0. 



Theorem 6.1 Lei < 6 < a be given. For uo G Q m ,n,{b 2 ,a 2 ] let C(u) be 
a positive ra x n matrix defined above and A(uS) G A m+n be given by \6. II) . 
(C = C(lo) ). Assume that Xij, 1 < i < j < (m + n) are independent Gaussian 
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satisfying U.S\) . Let Ya £ An be the random skew symmetric matrix given by 
(TO)) . Then for any t > 



limsup sup Pr( — — | logdet(v / t/j V +y A )-log $(£, !£„,,«,«) | > 5) = 
n^oo m<ti,wen miBj62o2] m + n 



(6.7) 

Equivalently 

limsup sup — — (log$(t,ir TOi „ )£l ,) -log$(t,K min;U; )) = 0. (6.8) 

Proof. Our proof follows the arguments in [9], and we point out the 
modifications that one has to make. Let N = m + n. Since 1 < m < n we 
have that ^ < ± < ±. ffl~2j) with = 1 implies: 



limsup sup Pr( — ^— | logdet(v / t/Ar+Y A )-log$(t, K m>n ^)\ > 8) = 

— 3 m,n, [o^ ,a^J 

(6.9) 

Thus it is enough to show equality (16.81) . 

Denote by Xa the random hermitian matrix Xa '■= -^=iYa- For e > 
define 

N 



det € (VtI N + Y N ) ■=Y[yt + vaax(\Xi{Y N )\, VNe) 2 , 
i=i 

ri N fl 

det e (^=I N - iX N ) := J] y - + max(\Xi(X N )\,e) 2 . 

Clearly, 

det e (Vt/iv + Vat) = det e (-^/ iV - iX N ). (6.10) 
Let fN,tA x ) '■= | lo g(^ + max(|x|,e) 2 ). Then 

|/iv,t, e U<-foriV>4- 
In what follows we assume that iV > -4. Observe next that 

-J- logdet e (-^=/AT - iXat) = tr N f NA e(X A ). 
t v v iV 



Combine the concentration inequality (13.11) with (16.101) to obtain 

1 r- r- N 2 r 2 e 2 

Pr(| — (logdet e (v / t/7v+^v)-E logdet e (Vt/jv+>Jv))| > r) < 2eSsr- (6.11) 
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Let 



(logiV) 2 



Note that for a fixed t one has iV > 4- for N » 1. Hence 



limsupPr(i-| logdet £iV (v / t/iv + Ijv) - E logdet^v^/jv + 1^)1 > 6) = 

for any 5 > 0. As in [9j Prf. of Lemma 2.1] use (16. lip and Lemma 13.21 to 
deduce that 

lim — {logE det eff (Vil N + Y N )-E \ogdet 6N (Vil N + Y N )) = 0, 

A^— >oo iV 

which is equivalent to 

lim l(logEdet ejv (^I iV -iX i v)-E logdet £N (^Ll N -iX N )) = Q. (6.13) 



It is left to show that under the assumption of the theorem 

lim ^(logE det eN {VtI N + Y N ) - logE det{Vil N + Y N )) = 0. (6.14) 



Clearly, the above claim is equivalent to 

lim ^-MogE det 6N (^Ll N - xX N ) - logE det(^Ll N - iX N )) = 0. (6.15) 



To prove the above equality we use the results of [9]. First observe that 
has at least n — m eigenvalues which are equal to zero, while the other 2m 
eigenvalues are ±Xi(X N ), . . . , ±X m (X N ). Furthermore Xi(X N ) 2 , . . . , A^(X/v) 
are the m eigenvalues of jjUcUc , denoted in [9] as Z(A n<m ). Clearly 

ft ft m t 

det e (^=I N - iX N ) = (^=)^ m II(^ + ™x(HXN) 2 ,ef) > 



det(^Ll N - iX N ) = (^) n - m fl(^ + \(X N ) 2 ). (6.16) 



i=i 

Hence for e < 1 

< ±(\ gdet e (^=I N - iX N ) - logdet(^J iV - iX N )) 
v N v N 



1 , ¥ + e2 1 v- , e2 

iv ^ log xJx^y' 
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[91 (3.2)] is equivalent to 

limsup sup E — — lo SY7^~ — ^2 = °- 

n-oo m<n,u>en [b2 , a 2j m + n x(x , \2 <f 2 M^m+n) 

Hence 

lim i-(E logdet eN (^Ll N -xX N )-E l gdet(^Ll N - iX N )) = 0. (6.17) 



Combine (16. 16ft with Jensen's inequality to deduce 

E \ogdet(^Ll N -iX N ) < logE det(^Ll N -iX N ) < logE det e (-^=I N -iX N ) 
Hence 

limsup -J=(logE det tN (^Ll N - iX N ) - E logdet(-^L/jv - iX N j) > 

n^oo V-/V V N 



limsup ^(logE det eAr (-^/iv - iX N ) - logE det(-^Lijv - iX N )) > 0. 
Use (I6U3|) and (IBTTTD to deduce (16TT5|) . □ 
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